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The nontrivial evolution of Wannier functions (WF) for the occupied bands is a good starting 
point to understand topological insulator. By modifying the definition of WFs from the eigenstates 
of the projected position operator to those of the projected modular position operator, we are able 
to extend the usage of WFs to Weyl metal where the WFs in the old definition fails because of the 
lack of band gap at the Fermi energy. This extension helps us to universally understand topological 
insulator and topological semi-metal in a same framework. Another advantage of using the modular 
position operators in the definition is that the higher dimensional WFs for the occupied bands can 
be easily obtained. We show one of their applications by schematically explaining why the winding 
numbers usd = h> 2 D for the 3D topological insulators of Dill class presented in Phys. Rev. Lett. 
114, 016801(2015). 
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I. INTRODUCTION 

The discovery of quantum Hall effect [l|, for the first 
time, sheds light on the topological properties of elec¬ 
tronic bands hidden behind the 2-dimensional (2D) Fermi 
gas in magnetic field. Thouless et al. (referred as TKNN) 
calculated the Hall conductance by Kubo formula [2| and 
expressed it as an integral in the whole first Brillouin zone 
(FBZ). Then Simon et al. recognized the integral as the 
first Chern class of the U(l) principal fiber bundle on 
the torus of the FBZ Besides this route, Laughlin 

raised another way to understand quantum Hall effect by 
a gauge invariant argument (referred as the Laughlin ar¬ 
gument on a cylinder) [6| . These two ways are equivalent 
but the latter one is more profound in explaining the ro¬ 
bustness of the Hall conductance against local disorders 
because it is based on the global gauge arguments of the 
2D electron system. 

For t(mological insulators (TI), e. ^ 2D spin Hall 
system Q and 3D strong or weak TI |8l-[lQ|. there are 
also two routes, parallel to those for quantum Hall ef¬ 
fect, to understand their topological properties. One is 
the topological invariant, likes the Hall conductance in 
the TKNN theorem, that can be expressed as an integral 
in a portion of the FBZ [ilil. The other route, simi¬ 
lar to the Laughlin argument, is based on the evolution 
of the Wannier centers (WC) of the localized Wannier 
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functions (WF) (along one direction) as a function of the 
momentum (in the perpendicular directions). These two 
routes are also equivalent. For instance, in Ref. Q, the 
Z 2 topological invariant in the 2D spin Hall system is 
defined by the switch partners of the WFs, which can be 
schematically illustrated by the evolution of WCs in the 
X direction as a function of momentum in the y direction, 
ky. In Ref. [il-[i3> the topological invariant for 2D and 
3D TIs are also associated with the evolution of the WCs 
as a function of ky (and kz) as well. Beside the schematic 
illustration of the nontrivial topological bands, WFs are 
also a good one-particle basis to study the many-body 
wave functions of the fractional topological insulator by 
mimicking the combination of these WFs as that of the 
Landau orbitals playing roles in the Laughlin state in 
the fractional quantum Hall effect [H, [TqI , or introduc¬ 
ing pump process to axion coupling directly. [ 2 ^ 

Besides TI, another group of topological materials, 
called topological semi-metal (TM), such as Weyl metal, 
has attracted great interest for their exotic surface states 
and bulk behaviors in recent years [2ll - [^ . Weyl metal 
does not have topological invariant defined in the 3D FBZ 
because bands touch at the Weyl points so that there 
is no band gap separating the occupied and the empty 
bands. One usually needs to split the 3D FBZ into a se¬ 
ries of 2D slides and the topological invariant, like Chern 
numbers, can only be defined on these 2D slides that do 
not cross the Weyl points. It is also impossible to directly 
adopt the evolution of WCs to understand TM because 
the discussion of WCs has a prerequisite condition: the 
WCs of WFs must be meaningful. But as discussed by 
Kivelso et ah, the WFs became extended as long as the 
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band gap at the Fermi energy became zero [23-1^. So 
the traditional WCs of WFs become meaningless in Weyl 
metals. 

In this paper, we adopt the modular operators, e.g. 
X = Xjnod + where x is the position operator in 

the X direction and Xmod is the modular position oper¬ 
ator, which can be expressed as Xmod = xmodL in the 
representation of x [^. Here L is the length of a unit 
cell. We change the definition of WFs, |WF), from the 
eigenstates of PxP to those of PxmodP^ where P is the 
project operator on the occupied states. The WCs of 
these WFs measured in each unit cell are the eigenval¬ 
ues of the projected modular position operator. We find 
that the standard deviation of projected modular posi¬ 
tion operator, (WF|P;r^^^P|WF) is zero for every WF 
IWF) even when the gap between the occupied bands and 
the empty bands goes to zero. As a result, the WCs of 
these newly defined WFs in each unit cell are still mean¬ 
ingful in both TM and TI. By this way, we can univer¬ 
sally illustrate the topological properties of TI and TM 
by evaluating the evolution of WCs. 

Although the physical model may be in 2D and 3D, 
the WFs defined above are 1-dimensional (ID) as they 
are the eigenstates of the operator Xmod and the mo¬ 
menta ky and kz are only taken as parameters. We find 
that such WCs of the ID WFs are equivalent to the phase 
accumulated by the Bloch states along the Wilson loop 
in the k^ direction, which has been introduced by Yu et 
al. and Soluyanov et al. individually One may 

regard the extension of the ID WFs as a new representa¬ 
tion of their conclusions, but we should emphasize that 
beside the application of the ID WFs to TM, from the 
present starting point the 2D or higher dimensional WFs 
can be easily obtained. Previously, such higher dimen¬ 
sional WFs are seldom used in TI due to the reasons such 
as [PxP, P^P] 7 ^ 0 that causes the absence of 2D WFs 
as common eigenstates of PxP and PyP. As a contrast, 
we will show that such obstacle is removed in our con¬ 
sideration because [PxmodPi PVmodP] = 0. Here y and 
ymod are the position and modular position operators in 
the y direction, respectively. 

This paper is organized as the following: In section im 
we show that the states in each band are also the eigen¬ 
states of the modular position operators with a certain 
degenerated eigenvalue Xmod- We then show that 3D TI, 
as well as 3D TM, can be universally illustrated in the 
same framework with {xmodi^y^kz). In section Hill we ex¬ 
tend our WFs from ID to higher dimensions. We show 
their applications by illustrating in a plausible way that 
why the winding numbers tj 2 D = in two models of 
the Dili class. In section IIYI we present the conclusions. 


II. ID MODULAR WANNIER FUNCTIONS IN 
MULTI-BAND SYSTEMS. 

Traditionally, the ID WFs for the occupied bands are 
defined as the eigenstates of the projected position oper¬ 


ator, 

X = PxP, (1) 

where x is the position operator in the x direction and 
is the project operator on the oc¬ 
cupied states in the Hilbert space. Here, ID refers 
to the direction of position operator x, while the mo¬ 
menta along the other dimensions in the 3D space of the 
model are taken as parameters. Without loss of gener¬ 
ality, we take the Fermi energy as energy zero. These 
WFs and the corresponding WCs can be denoted as jx^) 
and Xi, respectively. The second moment of the WF 
P = {xi\{x — Xi)‘^\xi) ^ where Eg is the band gap at 
the Fermi energy, will diverge as long as ^0. This is 
why there is no discussion made directly from the point 
of view of WFs for TMs. 

Next we show that this obstacle can be removed when 
the modular position operator is introduced. We will 
first introduce the modular variances in a ID continuum 
system and then applying them to ID discrete multi-band 
models. 

For a ID infinite continuum system in the x direction, 
one can introduce the modular position operator Xmod 
and the modular wave-vector operator kmod = Pmod/^ 
as 

X — NxL T XfYiQd (2) 

^ .. 27r 

k — P ^mod’) (3) 

where Pmod is the modular momentum operator in the x 
direction and Xmod (kmod) can only take the eigenvalues 
in range of [0,P) ([0,^) ) in the x (k) representation 

and Nx i^k) can only take integer values [28|. Here L is 
an arbitrary length. 

The system can be translated L in the x direction by 
applying a translation operator Simi¬ 

larly, to translate the system ^ in the momentum space, 
the translation operator needed. It 

is easy to find that ^g ^ikmodL 

— ^ ig a one-to-one mapping for kmod (xmod)^ we 
conclude that [xmod^ kmod] = 0 [28|. This means that the 
modular variances commute although the original vari¬ 
ances have [x^k] = 

Such modular variances can be naturally extended to 
the discrete periodic multi-band models in solid. We 
take L as the length of a unit cell so that kmod take 
the value in the FBZ and G = ^ is the reciprocal 
vector. We denote the ith orbital state in the nth 
unit cell as \n)i and the ath band state as \k)a, with 
\k)a = En=i E^i Uaj{k)\n)j, where N is the 
number of the unit cells and U{k) is a unitary matrix. 
There are totally m orbital states in each unit cell so 
that there are m bands in the FBZ. The wave-vector op¬ 
erator k and the position operator x are defined in the 


real space by 
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= \n- l)i (4) 

x\n)i = {nL + Xi)\n)i, ( 5 ) 

where Xi G [0, L). As and 

in the reciprocal space, by applying the above equations 

to \k)a, we have 

^ikrnodL^J^^^^ = e^^^\k)a (6) 

^ (7) 

In the last equation, we have used the property that 
in the reciprocal space, the states \k)a and \k + G)a 
must represent the same quantum state with only a pos¬ 
sible phase difference Oa G [0,27r). Here 9 a is inde¬ 
pendent of k because these transformations for different 
\k) will map to the same Wilson loop with only start¬ 
ing point difference. Like that in the continuum model, 
^Q^kmodL^^iXmodGj = Q go that krnod and Xmod are diago¬ 
nalized simultaneously in the reciprocal space for each 
band. These Wfold degenerated eigenvalues of Xmod 
for each band correspond to N WFs with their WCs at 
Xmod = 9oc/G in each unit cell. 

As Xmod is diagonal for each band, one can apply the 
occupied project operator P from the left and the right 
sides straightforwardly and the second moment of the 
eigenstates of the projected modular position operator, 
{xi\Px^^^j^P\xi) - {xi\PxmodP\xiP is fixed at zero. We 
want to emphasize that in the above derivations, we only 
suppose that each band is either fully occupied or empty 
but do not require a finite band gap at the Fermi energy. 

In the case of multi-bands being occupied, it will be 
hard to distinguish the band index a in the presence 
of band crossing. This difficulty can be overridden by 
studying the matrix This operator repre¬ 

sents how much phase is accumulated for the states \k)a 
by mapping to \k P G) a = |^)a- One can divide this 
mapping into Nk segments, with 

n/c G (1, 2, • • • , XP). Then the phases XmodG can be ob¬ 
tained by finding the phases of the eigenstates of the 
product matrix The above mapping is equiv¬ 

alent to the product of the non-abelian Berry connection 
along the “Wilson loop” raised in Ref. [l^, [iQ • l^^ie 
phases XmocjG = 9a are equal to the phases denoted by 
9^ in Ref. jl^ . 

Now we apply the above results to a 3D Weyl metal 
with the effective Hamiltonian 

Hw = [Mo/ 3i{cos{kx) P cos{ky)) P P2Cos{kz)]az 

+ sm{kPaa: + sm{ky)ay, ( 8 ) 

where is the Pauli matrix along x(yz) direction. 

Here we fiFft = 0.2 and ^2 = 0.4. 

The model has two Weyl points at {k^^ky^kz) = 
(0,0, zbl.29) when Mq = —0.5. Like the applications 
in 3D TI [l3, we can take the two wave-vectors, e.g. 
kx and kz, as parameters and find the eigenvalues of 



FIG. 1: 9y as Si function of kx and kz- 9y acquires a 27r (0) 
variation when {kx,kz) evolves along C 2 (ci). Here Nk = 20. 

PymodP = ^y/G along the rest y direction. In Fig. [H 
we plot the function 9y{kx, kz) in the range [—tt, tt]. 9y is 
27r folded along the line connecting the two Weyl points 
(projected on the kx — kz plane) which is the projection 
of the Dirac string on the kx — kz plane. Imagine that we 
are varying the wave-vectors along the typical paths Gi 
and G 2 shown in Fig. [TJ The WCs of WFs will come back 
to their original positions when the varying is along Gi. 
But the WCs change one unit cell when the varying is 
along G 2 ‘ By applying Laughlin’s cylinder argument, we 
can conclude immediately that there should be bound¬ 
ary states that cross the zero energy on the xz surface 
to compensate the movements of WFs to make the situa¬ 
tion totally unchanged when the wave-vector comes back 
along G 2 on the 2D FBZ torus. But such zero energy sur¬ 
face state does not exist when the varying is along Gi. 
So these zero energy boundary states form the unclosed 
Fermi arc on the x — z surface of the 3D Weyl metal. 

We find a similar plot for 9x{ky^kz)^ which indicates 
that Fermi arc can also appear at the y — z surface. But 
we find a trivial 9z{kx^ky)^ keeping at 0 for any kx and 
ky^ which indicates that there is no Fermi arc at the xy 
surface of the Weyl metal. 

The model can have 6 Weyl points at (0, tt, ±0.73), 
(tt, 0, ±0.73) and (0, 0, ±1.82) when Mq = —0.3. We plot 
the function 9y{kx,kz) in Fig. [2l As that in Fig. [H 
the plot illustrates the topological properties of the Weyl 
metal explicitly. The evolution of 9x{ky,kz) is also simi¬ 
lar to 9y(kx^kz) and 9z{kx^ky) is trivial as well. 

By the two illustrative models, we show the power of 
the WFs defined by modular position operator. For the 
first time, we can understand TM by considering the 3D 
FBZ as a whole, otherwise one usually needs to employ a 
series of 2D planes or 2D cylinders in the FBZ to discuss 
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FIG. 2: 6y as a function of kx and kz for the model with 
totally 6 Weyl points. 


topological properties of TM. 


III. EXTENDED WANNIER FUNCTIONS TO 
HIGHER DIMENSIONS 

In this section, we show that with the help of the mod¬ 
ular position operators, we can easily find WFs in higher 
dimensions. We will show the application of 2D WFs in 
two 3D illustrative lattice models in the Dill class. 

By applying Eq. [71 we have 


(9) 

y/ai (10) 


piXmodGxlU h \ — \ h h \ 

^iy,no,Gy\k^^ky)c. = 

where \kx^ky)a denotes the eigenstate in the ath band 
with momenta kx and ky in the 2D FBZ. Here, 

()a^^\ky(x)) is a function of ky^^) so that Xmod and ymod 
(the eigenvalues of the modular operators) are not uni¬ 
form in each unit cell. This is the key difference between 
WFs in higher dimensions and those in ID. 

As the modular position operators, Xmod and ymodi are 
diagonal for each band, we have [PxmodP^ PymodP] = 0. 
But for the position operators, one has [PxP, PyP] ^ 0 
[ 2 ^. This is why after defining the WFs by the modular 
position operators, we can extend WFs to higher dimen¬ 
sions easily. 

We next show how to find the WCs of the WFs in 
2D. The method is straightforwardly to be applied in 
further higher dimensions. We take a square 2D FBZ as 
an instance. We discretize the FBZ into N‘^ elements by 
N rows and N columns and denote the state in each cell 


as 


nx,ny)a = I - ^ + rixSkx,-^ + ny6ky)a, where 
Gx(y) is the reciprocal vector and Skx(y) = Gx(y)/N. We 
calculate the operators 


XY~^ = p 

XY~ = P ^ 


( 11 ) 

( 12 ) 


in the discrete FBZ. Most elements of XY^ are zero ex¬ 
cept X^nx,ny-,nx + l,ny + l and + — 1 

for XY~. Furthermore, XY^ can be block diagonalized 


into N independent N x N matrix, with each matrix 
mimics the Xp{ky) matrix in Eq. 6 in Ref. [TgI, 


/ 


XY^{n) 


XY, 


l,n;2,n±l 

0 


0 


AW 


2,n±l;3,n±2 




n± V;!,?! 


0 


0 


(13) 

The phases of the eigenvalues of the product matrix 
^y'!n;2,n±l^W±l;3,n±2 ■ • ■ ^y^,n±JV;l,n gi^C the N- 
fold degenerated values XmodGx^ymodGy (mod 27r). Like 
the ID case, these A-fold degenerated modular positions 
correspond to WCs localized at different unit cells. But 
as there are N independent AW^(n), we will have totally 
N kinds of A-fold degenerated values XmodGx i ymodGy 
(mod 27r). So we can conclude that the 2D WCs defined 
by these modular position operators are not uniform in 
each unit cell in the 2D lattice. This is consistent with the 
previous argument that the 2D WEs with uniform WCs 
do not exist in 2D topological nontrivial lattice [25|. 

Now we apply the above method to find 2D WEs for 
two 3D lattice models. We will take the momentum, 
e.g. kx^ as a parameter and study how the 2D WCs 
evolve as kx is varying 27r in the EBZ. Such illustration 
provides more information than the evolution of ID WCs 
as varying, e.g. both kx and ky. 

We start from a 3D illustrative model in the Dill class. 


- Wik) e{k)l) ’ 


(14) 


with 0(^) = i{d'd)(7y^ where the unit matrix / and Pauli 
matrix a = (o^cr^, (Tz) are all 2 x 2 matrix [^ . 

As in Ref. [^, we denote the Hamiltonian as iLi when 
e(A?) = cos{kx — kz)-\-cos{kx)-\-2cos{ky — kz)-\-cos{ky)—ii 
and {dx^ dy^dz) = A([sin(A:a^ — kz) Y sm{kx) Y 2 sm{ky — 
kz)—sm{ky)], 2 sin(/cy), 2 sm{kxYky)). The winding num¬ 
ber Usd for this model can take nontrivial value ±1 [29|. 
Here we take A = 5 and fi = 2.5 so that usd = 1 in our 
calculation. Einch et. al argued in Ref. that when 
a breaking symmetry term is added to the Hamiltonian 
on the surface of this TI, e.g. on the xy surface, to open 
up a gap for the surface states, one is able to obtain the 
winding number for the gapped surface states, U 2 Di and 
found that usD = 1^20 > 

In Eig. [3] (a) and (b), we illustrate how O^Pn) = 

GiUmod^Zmod) evolve with where G = = Gy = ^ 

is the reciprocal vector and ymodpmod) is the eigenvalues 
of the modular position operators ymodpmod)- Here, for 
the sake of clarity, we divide the whole EBZ into 4x4 
blocks so that there are 4 kinds of 0^^{n), which are iden¬ 
tified by the index n = 0,1,2,3 in the figure. Dividing 
the EBZ into smaller parts will not alter the key features 
of the figures. As the figure shows, most of the modu¬ 
lar positions experience a trivial evolution with kx- But 
there is one modular position (one for ymod + ^mod and 
one for ymod — ^mod respectively) that undergoes a change 
of unit cell a as varying a reciprocal vector for kx. 
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FIG. 3: O^z — G{ymod^Zmod) as functions of kx in (a) and (b). 
The 2D FEZ in ky and kz plane is divided in to 4 x 4 blocks, 
n = 0,1, 2, 3 are the 4 kinds of O^z- Here the lattice constant 
a is taking as the unit, (c) shows the 3D topological insulator 
with breaking symmetry term added on the xy surface. We 
show how the WCs of WFs evolve with kx in the 2D real space 
in (d), with red and green dots representing the two WCs in 
each unit cell. The arrows indicate the variations of WCs by 
changing kx from 0 to 2'k ja. Most of them undergo a trivial 
evolution except those along a column and a row. 


Such behavior of evolution reflects the bulk topologi¬ 
cal properties of the 3D model which can not be modified 
by adding a symmetry breaking term on the xy surface. 
We explicitly illustrate how these WCs evolves with kx 
in the 2D real plane yz in Fig. [3] (d). Points in each 
unit cell are representing the WCs of the WFs. The 
exact order of columns or rows are not important be¬ 
cause a gauge transformation can exchange them. This 
figure shows that there is one WF that being pumped 
one half lattice constant from left to right (and another 
WF that is pumped one half lattice constant inversely) 
along y axis as varying kx circularly from 0 to 2^ ja. To 
compensate such evolution for the gauge invariant, there 
should be topological protected edge states at the edge 


of either gapped or gapless xy surface. Because of the 
bulk-boundary correspondence for TI, each xy surface 
should contribute 1/2 winding, which takes the respon¬ 
sibility of the half lattice movement of the WFs during 
the kx pump. As a result, which counts the total 
winding of the two xy surfaces, is totally 1. This is what 
= ^ 21 ) means for this model. 



FIG. 4: O^z = G{ymod + Zmod) as functions of kx in (a) and 
(b). The 2D FEZ in ky and kz plane is divided in to 4 x 4 
blocks, n = 0,1, 2, 3 are the 4 kinds of 6^z- 

We then do a similar discussion on another model H 2 
with e{k) = cos{kx) + cos{kx — kz) — yu and {dx^dy, dz) = 
A{sm{kx) + sm{kx — kz),2sm{ky),2sm{kx — kz)). We 
take A = 2 and /i = 1 in the calculation. The model is 
in the topological phase with = 2. 

In Fig. m we show how G{ymod + ^mod) and G{ymod - 
Zmod) evolve with kx. There are two modular positions 
experience nontrivial evolution along ^ -h ^ or ^ — 2 : direc¬ 
tion. This makes us infer that there are two WCs being 
pumped half a unit cell from left to right and another two 
WCs vice versa. Applying the similar discussion as we 
have done for the previous model we can conclude that 
U 2 d = 2. This also matches the equation = ^ 2 D for 
this model. 

We have also applied similar calculations for Hs and 
H 4 in the topological phase with 1 ^ 30 =^ and usjj = 4 
1^ . But we find that the evolution of 2D WCs are not 
simply extended as that from Hi to H 2 . Typically for 
i^ 4 , we find the evolution of 2D WCs are all trivial. One 
possible explanation is that the evolution of 2D WFs can 
not universally represent all TIs of the Dill class. But 
we propose that it is because another factor which has 
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not been discussed in Ref. [29|, the disorder effect. Our 
discussion is based on a gauge argument. So the conclu¬ 
sion naturally possesses the effect of disorder even though 
the disorder term is not explicitly engaged. For i74 with 
= 4, there is two pair of chiral edge states at the 
boundary of the 2D xy surface. Although the symmetry 
may still prevent the intra-pairing back scattering, it can 
not prevent the inter-pairing back scattering in the pres¬ 
ence of disorder. So the disorder can smear out the edge 
states near the zero energy and open up an effective gap. 
This is why we find the trivial 2D WCs evolution for 7^4 
and it will need further numerically investigation to be 
proved. 

IV. CONCLUSIONS 

We modify the definition of WFs for the occupied 
bands by introducing the modular position operators. 


One can calculate these ID WFs along the Wilson loop 
along one direction of the FEZ. The topological nontriv¬ 
ial band structure can be illustrated by the nontrivial 
evolution of WCs as varying the momentum in the other 
directions. We show that such illustration, usually only 
applied to TI, can also be applied to Weyl metal in which 
there is no band gap at the Fermi surface. Our illustra¬ 
tion well explains the nontrivial band structure of Weyl 
metal as well as the Fermi arc for the surface states. We 
also extend the ID WFs to higher dimensions. We show 
its application on explaining 1 ^ 30 = ^ 2 D in the two mod¬ 
els belonging to the Dill class. 
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